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The properties of the thermal Abehan color-magnetic monopoles in the maximally Abelian gauge 
are studied in the deconfinement phase of the lattice SU (2) gluodynamics. To check universality of 
the monopole properties we employ the tadpole improved Symanzik action. The simulated annealing 
algorithm combined with multiple gauge copies is applied for fixing the maximally Abelian gauge to 
avoid effects of Gribov copies. We compute the density, interaction parameters, thermal mass and 
chemical potential of the thermal Abelian monopoles in the temperature range between Tc and 3Tc. 
In comparison with earlier findings our results for these quantities are improved either with respect 
to effects of Gribov copies or with respect to lattice artifacts. 
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I. INTRODUCTION 

The signatures of the strong interactions in the 
quark-gluon matter were found in heavy ion collisions 
experiments [l]. There are proposals [1, Q suggest- 
ing that color-magnetic monopoles contribution can 
explain this rather unexpected property. These pro- 
posals inspired a number of publications devoted to 
the properties and possible roles of the monopoles in 
the quark-gluon phase P-[Toj . 

Lattice gauge theory suggests a direct way to 
study fluctuations contributing to the Euclidean space 
functional integral, in particular, the color-magnetic 
monopoles can be studied. In a number of pa- 
pers the evidence was found that the nonperturba- 
tive properties of the nonabelian gauge theories such 
as confinement, deconfining transition, chiral symme- 
try breaking, etc. are closely related to the Abelian 
monopoles defined in the maximally Abelian gauge 
(MAG) [nl, 113. This was called a monopole dom- 
inance [ij]. The drawback of this approach to the 
monopole studies is that the definition is based on 
the choice of Abelian gauge. There are various ar- 
guments supporting the statement that the Abelian 
monopoles found in the MAG are important phys- 
ical fiuctuations surviving the cutoff removal: scal- 
ing of the monopole density at T = according to 



dimension 3 for infrared {percolating) cluster [14|; 
Abelian and monopole dominance for a number of in- 
frared physics observables (string tension (l3l - [lH , chi- 
ral condensate jlGj, hadron spectrum fT7l|); monopoles 
in the MAG are correlated with gauge invariant ob- 
jects - instantons and calorons [l^, [l^l- It has been 
recently argued that the MAG is a proper Abelian 



gauge to find gauge invariant monopoles since t'Hooft- 
Polyakov monopoles can be identified in this gauge 
by the Abelian flux, but this is not possible in other 
Abelian gauges ^2^. Most of these results were ob- 
tained for SU{2) gluodynamics but then confirmed 
for SU{3) theory and QCD [U, Listed above 

properties of Abelian monopoles survive the contin- 
uum limit and removal of the Gribov copy effects. It 
is worth noticing that removal of Gribov copy effects 
changes numerical values of monopole characteristics 
quite substantially [23|. 

In this paper we are studying thermal monopoles. 
It was shown in Ref. [3j that thermal monopoles 
in Minkowski space are associated with Euclidean 
monopole trajectories wrapped around the tempera- 
ture direction of the Euclidean volume. So the den- 
sity of the monopoles in the Minkowski space is given 
by the average of the absolute value of the monopole 
wrapping number. 

In 0] another approach to study thermal monopole 
properties in the quark-gluon plasma phase based on 
the molecular dynamics algorithm was suggested and 
implemented. The results for parameters of inter- 
monopole interaction were found in agreement with 
lattice results 0- 

First numerical investigations of the wrapping 
monopole trajectories were performed in SU(2) Yang- 
Mills theory at high temperatures in Refs. |24| 
and [2^. A more systematic study of the thermal 
monopoles was performed in Ref. ^S]- It was found 
in Q that the density of monopoles is independent 
of the lattice spacing, as it should be for a phys- 
ical quantity. The density-density spatial correla- 
tion functions were computed in It was shown 
that there is a repulsive (attractive) interaction for a 
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monopole-monopole (monopole-antimonopole) pairs, 
which at large distances might be described by a 
screened Coulomb potential with a screening length 
of the order of 0.1 fm. In Ref. it was proposed to 
associate the respective coupling constant with a mag- 
netic coupling am- In the paper Q trajectories which 
wrap more than one time around the time direction 
were investigated. It was shown that these trajectories 
contribute significantly to a total monopole density at 
T slightly above Tc- It was also demonstrated that 
Bose condensation of thermal monopoles, indicated 
by vanishing of the monopole chemical potential, hap- 
pens at temperature very close to Tc. However, the 
relaxation algorithm applied in [^] to fix the MAG is a 
source of the systematic errors due to effects of Gribov 
copies. It is known since long ago that these efi^ects 
are strong in the MAG and results for gauge nonin- 
variant observables can be substantially corrupted by 
inadequate gauge fixing ^J. For the density of mag- 
netic currents at zero temperature it might be as high 
as 20%. 

For nonzero temperature the effects of Gribov 
copies were not investigated until recently. In a recent 
paper flO*! this gap was partially closed. It was shown 
that indeed gauge fixing with SA algorithm and 10 
gauge copies per configuration gives rise to the den- 
sity of the thermal monopoles 20 to 30% lower (de- 
pending on the temperature) than values found in . 
Large systematic effects due to effects of Gribov copies 
found in Ref. jlOi] imply that results obtained in earlier 
papers 0-13 for the density and other monopole prop- 
erties can not be considered as quantitatively precise 
and need further independent verification . The quan- 
titatively precise determination of such parameters as 
monopole density, monopole coupling and others is 
necessary, in particular, to verify the conjecture [3] 
that the magnetic monopoles are weakly interacting 
(in comparison with electrically charged fluctuations) 
just above transition but become strongly interacting 
at high temperatures. In this paper we use the same 
gauge fixing procedure as in Refs. to avoid sys- 

tematic effects due to Gribov copies. 

The careful study of the finite volume and finite lat- 
tice spacing effects was made in [6|] . We fix our spatial 
lattice size to Ls — 48 which was shown in Ref. Q to 
be large enough to avoid finite volume effects. We 
check finite lattice spacing effects comparing results 
obtained on lattices with Nt — 4: and 6 at two temper- 
atures. Let us emphasize that our studies are compu- 
tationally much more demanding in comparison with 
studies undertaken in Refs. [1, since we produce 
10 Gribov copies per configuration to avoid Gribov 
copies effect. For this reason our check of the contin- 
uum limit is not as extensive as it was in Refs. 0, |^. 

The important contribution of this work to the ther- 
mal monopole studies is a check of universality. In 
studies of magnetic currents at zero temperature it 
was found ^4j that the density of the infrared mag- 



netic currents is different for different lattice actions 
with difference as large as 30%. The conclusion was 
made that the ultraviolet fluctuations contribute to 
the infrared density and this contribution has to be 
removed. Partial removal was made by the use of 
the improved action. In present paper we use the im- 
proved lattice action - tadpole improved Symanzik ac- 
tion and compare our results for the density and other 
quantities with results obtained with the Wilson ac- 
tion '6'-'8','To'|. We flnd that the universality holds for 
monopoles which do not form short range (ultraviolet) 
dipoles. 

We also want to point out that in this work we use 
more natural procedure of computing monopole cor- 
relators in comparison with papers [8|] and [l3]. It 
was mentioned in Ref. Q that monopole trajectories 
had a lot of small loops attached to them which were 
UV noise. Presence of such loops do not allow to 
determined monopole spatial coordinates unambigu- 
ously for all time slices. This problem was bypassed 
in [l3| by using only one time slice. We remove the 
small loops attached to thermal monopole trajectories 
and thus we are able to determine the monopole co- 
ordinates in every time slice unambiguously. Then we 
use all time slices to compute the correlators what al- 
lows us to decrease the statistical errors substantially. 



II. SIMULATION DETAILS 

We studied the SU{2) lattice gauge theory with the 
tadpole improved Symanzik action: 



S — Pimpr ^ ^ 



pi 



trnpr 



(1) 



where Spi and Srt denote plaquette and 1x2 rectan- 
gular loop terms in the action: 



Spl.rt = 7;Tr{l - Upl^rt), 



(2) 



parameter uq is the input tadpole improvement factor 
taken here equal to the fourth root of the average pla- 
quette P = {\Upi). We use the same code to generate 
conflgurations of the lattice gauge field as was used 
in Ref. [ijj. Our calculations were performed on the 
asymmetric lattices with lattice volume V — LtL^, 
where Lt^s is the number of sites in the time (space) 
direction. The temperature T is given by: 



T 



1 

aLt 



(3) 



where a is the lattice spacing. To determine the val- 
ues of uq we used results of Ref. ']14i\ either directly or 
to make interpolation to necessary values of l3. The 
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critical value of the coupling constant for Lt = 6 is 
/3c = 3.248 j26|. For U = 6 the ratio T/T^ was ob- 
tained using data for the string tension from Ref. [T^ 
again either directly or via interpolation. For Lf = 4 
ratio T/Tc was taken to be equal to the ratio for 
Lt — 6 multiplied by factor 1.5. In Table |T] we provide 
the information about the gauge field ensembles and 
parameters used in our study. 

The MAG is fixed by finding an extremum of the 
gauge functional: 

Puid)^ i^E ^ Tr (u^^a^U^las) , (4) 

with respect to gauge transformations gx of the link 
variables Uxfi- 

Ux^ ^ U^^. = 9lUx^.9x+^ ; 9x e SU{2) . (5) 

We apply the simulated annealing (SA) algorithm 
which proved to be very efficient for this gauge |23l| as 
well as for other gauges such as center gauges [23| and 
Landau gauge To further decrease the Gribov 

copy effects we generated 10 Gribov copies per con- 
figuration starting every time gauge fixing procedure 
from a randomly selected gauge copy of the original 
Monte Carlo configuration. 
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TABLE I: Values of /?, parameter uq, lattice sizes, tem- 
perature, number of measurements used in this paper. 



III. MONOPOLE DENSITY 

The monopole current is defined on the links {x, fi}* 
of the dual lattice and take integer values j^(cc) = 
0,±1,±2. The monopole currents form closed loops 
combined into clusters. Wrapped clusters are closed 
through the lattice boundary. The wrapping number 
N^r € Z of a given cluster is defined by: 



jA[x)^cluster 

The density p of the thermal monopoles is defined 
as follows: 

_ ( '^clusters Wwr\ ) 

where the sum is taken over all wrapped clusters for 
a given configuration. 

Following Ref. [1] we distinguish two regions of tem- 
peratures: low temperature region T < 2Tc and high 
temperatures T > 2Tc. In Ref. [s^] it was proposed 
that at low temperatures the density of monopoles 
is almost insensitive to temperature thus indicating 
that the monopoles form a dense liquid while at high 
temperature the monopole density has a power depen- 
dence on temperature. These statements were based 
on the results obtained in Refs [l^l and [1^. The 
behavior of the density at low temperature was not 
discussed in earlier papers d, H, • 

The results for a range of temperatures Tc < T < 
1.5Tc obtained on lattices with Lt = 6 are presented 
In FIG. [TJ It can be seen that the thermal monopole 
density monotonously increases as temperature grows. 
The density at T/Tc = 1.5 is 2.3 times as large as 
one at T/Tc = 1-03 what indicates that the monopole 
density is considerably sensitive to temperature. This 
fact allows us to say that previous results on this ob- 
servable did not reflect real situation and the conclu- 
sion made in Ref. 0] was incorrect. In FIG. [T] we also 
show the results from Ref. for comparison. The re- 
sults demonstrate similar dependence on temperature, 
however our results are lower at any temperature. It 
is due to Gribov copy effect. 

In FIG. Owe show our data for temperatures T > 
l.bTc- The data of Refs. Q and [l3| are also shown 
for comparison. As was concluded in Ref. [l3| effects 
of Gribov copies in results of Ref. Q are between 20% 
at T/Tc = 2 and almost 30% ior T/Tc 7. Our re- 
sults deviate from those of Ref. [10| by about 10%. 
Thus we observe violation of universality of the ther- 
mal monopole density at given lattice spacing. We 
need to check whether the continuum limit is reached 
in case of the Symanzik action for our lattice spacings. 
To answer to this question we compare results ob- 
tained on Lt = 4 and Lt = 6 lattices at at T = 1.5 Tc 
and T ^ 2 Tc. One can see from FIG. H that the 
change of the density with decreasing lattice spacing 
is small for both temperature values. Quantitatively 
it is less than 1% for T = 2 Tc and about 3% for 
T = 1.5 Tc- This allows us to state that in case of 
Symanzik action the results for the monopole density 
obtained on Lt = 4 lattices can be considered as be- 
ing close to the continuum limit. Similar conclusion 
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FIG. 1: The behavior of the monopole density(normahzed 
by Tc ) at low temperatures(blue empty circles). The re- 
sults from Ref. Q are presented for comparison(red empty 
triangles). 



about closeness to the continuum limit was made for 
the Wilson action in Ref. jl^l- Thus the violation of 
universality seen in FIG. [5] may persist to continuum 
limit. Note that at zero temperature the univers ality 
breaking effects were found to be much stronger [14l |. 
up to 30%. 

Different values for monopole density in case of 
Symanzik and Wilson actions can be explained by a 
fact that configurations have different number of short 
range dipoles at the same temperatures. By short 
range dipoles we imply paired monopole trajectories 
having opposite wrappings and separated by distance 
of order of one lattice spacing. Thus, if we calculate 
monopole density omitting small distance, it should 
bring densities closer to each other. We found indeed 
that the distance between two densities decreases after 
we remove dipoles of size a. 

We will come back to discussion of the universality 
for the thermal monopole density in the next section. 

The dimensional reduction suggests for the den- 
sity p the following temperature dependence at high 
enough temperature: 



piT) 



1/3 _ 



Cpg'iT)T 



(8) 



where the temperature dependent running coupling 
g^(T) is described at high temperature by the two- 
loop expression with the scale parameter At'. 



9-\T) = 



11 
12^ 



17 

4471-2 



(ln(21n(r/AT)) (9) 



We fit our data to function determined by equations 
dS]) and dH]). The good fit with x^/dof ^ 0.28 was 



FIG. 2: The dependence of the thermal monopole den- 
sity(normahzed by T"^) on temperature(circles). The line 
is a fit to ec|.((9|. The data from Ref. j^] (triangles) and 
from Ref. [10] (squares) are presented for comparison. The 
densities at 1.5Tc and 2Tc with Lf = 6 are labelled by the 
crosses. 



obtained for T > 2Tc. The values for fit parameters 
were Cp = 0.160(6), At/Tc = 0.144(3). Thus for high 
enough temperature the density p{T) is well described 
by the form which follows from dimensional reduction. 
The values of parameters Cp and differ from values 
obtained in [10| though the difference in results for the 
density is small. But note that we used fit over range 
of temperatures between 2Tc and 3Tc while in [lo| T > 
3Tc were used. The data for the thermal monopole 
density are presented in TAB. |Vl 



IV. MONOPOLE INTERACTION 

We computed two types of monopole density cor- 
relators g{r), for monopoles having the same charges 
(MM correlator) and for monopoles having opposite 
charges (AM correlator). The correlators are defined 
as follows: 

(pM(0)pA/(r)) , {pAiO)pA{r)) 

gMM(r) = 2 1 2 (10) 



pIi 



{pA{0)pM{r)) , (pAf(O)pA(O) 

9AM{r) = \ (11) 



PaPm 



PaPm 



where pm,a{Q) and pM.Aij") are local densities. It can 
be reexpressed as: 



^^^^ , . ^ J 1 , rfiVM(r) , ^ _1 1 / dNAjr) ^^ 

Pm 47rr2 dr pA 47rr^ dr ' 

(12) 
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FIG. 3: The correlation function QAAiir) for monopole- 
antimonopole case for three values of temperature: STc, 
1.5c and l.OTTc with = 4, 4 and 6 respectively. 



FIG. 4: The correlation function gMM{r) for rnonopole- 
monopole case for the same values of temperature as in 
Fig. |3l The lines are fits for these three cases. 



where NMir) is the number of monopoles, and simi- 
larly for gAAi- 

In our computations following Refs. we take 

dN(r) to be a number of monopoles in a spherical 
shell of finite thickness dr — a a,t a distance r from a 
reference particle, whereas Airr^dr is equal to a volume 
of this shell, i.e. number of lattice sites in it. 

Correlators gMM,AM{j') were calculated for nine 
temperatures in the range between l.OSTc and STc. 
Three AM and MM correlators are presented in FIG.|3] 
and FIG. |4] respectively. 

One can see that in AM case the correlators start 
below 1, reach maximal value equal to 1.3 — 1.4 at 
distance r between 0.5/T and 0.8/T and then decrease 
down to 1 at distance between 1.8/T and 2.5/T. In 
other words the repulsion at short distances changes 
to attraction at large distances. In MM case a picture 
is different, the correlators start at small value below 
1 and smoothly reach 1 from below at same distances 
as AM correlators do. In this case we can see only 
the repulsion. These results are in a good qualitative 
agreement with results obtained in Refs. (gI [l0|. 

It should be mentioned that we applied more cor- 
rect procedure to compute MM and AM correlators. 
In Refs. [To'l and ^ only one time slice was used to cal- 
culate the correlators since the monopole coordinates 
could not be determined for all time slices unambigu- 
ously because of small loops attached to monopole 
trajectory. We used a special algorithm which cut off 
all UV loops out of each monopoles trajectory. Then 
we could take into consideration all time slices to com- 
pute the correlators. This leads to more precise results 
for the correlators. For example, at T = 3Tc the sta- 
tistical errors decreased by factor 1.7 approximately. 



One can also see that correlators in FIG. [3] and 
FIG.|3]are temperature dependent. It is possible to fit 
the data with the temperature dependent potential: 

9MALAM{r) = e-^f"^)/^ (13) 

where V{r) can be well approximated by a screened 
potential: 

V(r) ^ ^e"""'^ (14) 
r 

at large distances. In eq. (|T4l) m^) is a screening mass 
and am is a magnetic coupling. 

To decide which data points can be included in the 
fit range we used the following method. We plot the 
dependence of rV{r) on rT in log scale for all tem- 
peratures (see FIG.[S]). The linear dependence should 
be valid at distances where eq. (ITil) is satisfied. Thus 
data at small distances which do not fall onto a line 
are to be discarded. Typically we discarded three data 
points. Data at large distances were discarded because 
of high statistical errors. Our results for the fit param- 
eters are presented in Table |TT1 Note that our values 
for /Ndof shown in Table [II] are in general lower 
than values of this quantities characterizing fit qual- 
ity obtained in Refs. (7l.[To|. 

In order to check finite lattice spacing effects we 
compared correlators and respective fit parameters on 
lattices with Lf = 4 and 6 at T l.bT^ and 2Tc. 
The comparison of these correlators for T = 2Tc is 
presented in FIG. [6] One can see that for gMM{r) the 
finite lattice spacing effects are small at all distances 
being maximal at distances rT w 1.2. For 5am(?') the 
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FIG. 5: The dependence of rV{r) on rT. Note log scale 
for X axis. The red line is the fit to eq. (|14l) . 
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TABLE II: Values of the magnetic coupling am and the 
screening mass mu obtained by fitting of MM correlators 
to eq. (|13|l . The double line separates the parameters ob- 
tained for Lt = 4(above the line) from those obtained for 
Lt = 6(below the line). 



good agreement is also observed at distances to the 
right of the distance corresponding to the maximum 
of the correlator while at short distances the finite 
lattice spacing effects are large. Results for T = 1.5Tc 
are similar. 

In FIG. [7] we compare the MM and AM correlators 
computed in this work and in Ref. [l^] at T '--^ 1.5Tc. 
It can be seen that the correlators are in good agree- 
ment with exception for AM correlator at small dis- 
tances. This implies that we observe universality of 
the correlators apart from distribution of the small 
dipoles: with Wilson action we observe more such 
dipoles than with improved action. Above we have 



FIG. 6: Comparison of MM and AM correlators com- 
puted at T = 2Tc on lattices with for Lt = 6 (empty 
blue squares) and Lt = 4 (filled red circles). 
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FIG. 7: Comparison of MM and AM correlators com- 
puted at r = 1.5Tc in this work(red filled circles) and 
Ref. [l3|(blue empty squares). 



concluded that the number of small dipoles is decreas- 
ing with decreasing lattice spacing. Thus, we may ex- 
pect that in the continuum limit the universality will 
restore at small distances as well. 

The dependence of the fit parameters, mjj and am-, 
on temperature, obtained for MM correlators, is pre- 
sented in FIG. [5] and FIG IHl respectively. One can 
see that the behavior of the fit parameters at T close 
to Tc is different from their behavior at high tem- 
perature. Close to Tc, in the range between l.OSTc 
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and 1.1 Tc we observe decreasing of both mjj and am- 
Such behavior was not observed before. This indi- 
cates that just above the transition the monopole in- 
teraction becomes weaker with increasing tempera- 
ture. Since near the phase transition the finite volume 
effects might be strong this observation should be ver- 
ified on larger lattices. At a bit higher temperature, 
T = 1.2Tc, both parameters jump to higher values and 
then their dependence on the temperature becomes 
different. While for m^) we observe slow decreasing, 
in agreement with results of Ref. (lOj .magnetic cou- 
pling am is slowly increasing. This is in qualitative 
agreement with Refs. 0, [l^. Coming to quantitative 
comparison we find that our values for mo are almost 
within error bars although systematically lower than 
values reported in Ref. [l3|. The values for am pre- 
sented in FIG [9] are substantially lower than the val- 
ues obtained in Ref. In the temperature range 
T > 2Tc our results for am are also lower than val- 
ues obtained in Ref. Q . The data in FIG [S] indicate 
that am approaches its maximum value of about 3 at 
rather small temperature T ^ 2.5Tc. But for temper- 
atures close to Tc our values of am are higher than 
values presented in [7|. 

The difference in the behavior of the correlator pa- 
rameters can be explained at least partially by the fact 
that in this work we eliminated all small loops out of 
each wrapping cluster. This gave us an opportunity 
to determine the monopole location in a given time 
slice unambiguously and thus, to use all time slices 
for correlators computation. This procedure was used 
in studies of the thermal monopoles for the first time. 

We now can compute the plasma parameter F which 
is defined as follows: 



1/3 



(15) 



F is equal to ratio of the system potential energy to 
its kinetic energy. If F << 1 the system is a weakly 
coupled plasma; if F > 1, it is a strongly coupled 
plasma. For 1 < F < Fc 100, the system is in a liq- 
uid state. The dependence of F on temperature is pre- 
sented in FIG. 1101 F is roughly proportional to a since 
p^/^/T varies slowly with temperature (see FIG. [1] 
and [2]). Thus at small temperature we observe in 
FIG. [To] the nonmonotonic behavior we saw in FIG[9l 
At temperatures above l.STc F can be well approx- 
imated by a constant. We cannot exclude slight in- 
crease or slight decrease of F for higher temperatures, 
though. The independence of F on temperature was 
predicted in Ref. [7]. But quantitatively our result is 
rather different: in Ref. 7] F was found substantially 
higher and approaching a constant value about 5 at 
temperatures above 4Tc. Despite this quantitative dif- 
ference we confirm result of Ref. 7] that the thermal 
monopoles are in a liquid state at all temperatures. 
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FIG. 8: The dependence of the screening mass mo on 
temperature. 
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FIG. 9: The dependence of the magnetic couphng Om on 
temperature. 



V. MONOPOLE CONDENSATION 

In this section we consider thermal monopole trajec- 
tories which wrap more than one time in a time direc- 
tion. It was proposed in Q that these trajectories can 
serve as an indicator of Bose-Einstein condensation of 
the thermal monopoles when the phase transition is 
approached from above. The main idea of this pro- 
posal is the following: a trajectory wrapping k times in 
a time direction is associated with a set of k monopoles 
permutated cyclically. Having such an interpretation 
one can assess a density of these trajectories assuming 
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FIG. 10: The dependence of the Coulomb plasma param- 
eter on temperature. 



that they form a system of non-relativistic noninter- 
acting bosons. According to Ref. Q this density can 
be written as follows: 



-fik 



Pk 



(16) 



where is a number of wrappings, ft = -~^jl/T is a 
chemical potential, and A is the Dc Broglie thermal 
wavelength. The condensation temperature is deter- 
mined by the vanishing of the chemical potential. 

To take into account interactions between 
monopoles it was suggested in [8|] to modify eq. (HI 
to 



Pk 
T3 



(17) 



with a free parameter a. The condensation of 
monopoles still should be signalled by vanishing of 
effective chemical potential p. 

One can see from FIG. [TT] that our values for 
^ are systematically lower than values presented in 
Ref. This is in a qualitative agreement with the 
fact reported above that our total density p is sub- 
stantially lower than total density found in Ref. 
As in the case of the total density this difference is to 
be explained mostly by large systematic effects due to 
Gribov copies in results of Ref. . 

We fitted our data for ^ to eq. p7)) for tempera- 
tures close to Tc- Since our data do not allow us to 
keep free all parameters we made fits for 4 fixed values 
of a. The results of the fit are presented in Table Hill 
One can see that is decreasing with increasing a 
for three temperature values out of four. 



FIG. 11: Normahzed densities for trajectories wrapping 
more than once in a time direction (empty symbols). Note 
a log scale for Y-axes. For comparison we show here the 
data from Ref. [|| (filled symbols). 





l.OSTc 


1.07rc 


l.lTc 
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/i(a = 2) 


0.50(9) 
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1.6(1) 


2.0(3) 


M(a = 2.5) 


0.37(8) 


0.82(4) 


1.4(1) 


1.8(2) 


/i(a = 3) 


0.26(7) 


0.65(5) 


1.2(1) 


1.6(2) 


X\a = 0) 


2.09 


0.46 


2.27 


1.99 


X'(« = 2) 


1.35 


0.27 


1.84 


1.54 


x'(« = 2.5) 


1.22 


0.41 


1.73 


1.42 


X'(a = 3) 


1.14 


0.64 


1.61 


1.30 



TABLE III: Values of fit parameters p obtained by fitting 
data for ^ to eq. (|17|) for 4 values of parameter a. Values 
of are also shown. 



We also computed the monopole thermal mass using 
the mean squared monopole fluctuation Ar^ . The 
relation between Ar^ and the mass of a nonrelativistic 
free particle is as follows Q: 



1 



2TAr 



(18) 



Ar'^ can be computed on a lattice in the following way 



^Ar^ 



1 ^ 

1=1 



(19) 



where L - is a total trajectory lengths, df is a squared 
spatial distance between the monopole coordinate at 
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FIG. 12: The dependence chemical potential on temper- 
ature(blue filled triangles). The data from ref. @] is pre- 
sented for comparison (red empty triangles). 
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FIG. 13: The dimensionless ratio m/T as a function of 
temperature obtained in our work {Lt = 4 - red filled 
squares, Lj = 6 - red filled circles), in Ref. [lol| (empty 
triangles) and in Ref. Q, a = 0.047 (blue empty circles). 



i = and its current coordinate after i steps along the 
monopole trajectory. 

As in Refs. [13] we computed Ar^ for trajecto- 
ries with one wrapping. Furthermore, we have checked 
the influence of the small loops on the value of the 
monopole mass determined via eq. (jl8p . 

The comparison of the thermal monopole mass ob- 
tained for the Wilson action [l^] with the values 
obtained in this paper is presented in FIG. [13] It is 
seen that the values obtained for the Symanzik action 




FIG. 14: The dependence the average number of loops 
per one wrapped cluster on temperature. 



demonstrates the same dependence on temperature as 
for the Wilson case, but are higher at all temperatures 
than results of [I^ and lower than results of Q. 

The monopole mass computed after removal of the 
loops increases by a factor of 6% at high tempera- 
ture in comparison with the mass obtained in case 
when all loops are untouched. But the difference be- 
tween two masses is getting larger as temperature ap- 
proaches to Tc reaching 16% at l.OSTc- Such behavior 
of the thermal monopole mass is expectable as with 
decreasing temperature the number of loops attached 
to a wrapped cluster increases. This can be seen from 
FIG. [HI where the temperature dependence of num- 
ber of loops per one wrapped cluster is presented (rip). 
When temperature decreases from 3Tc to l.OBTc, dif- 
ference in rip is one order. 

It was found in 10] that the thermal monopole mass 
can be fitted by the following function: 



T A 



(20) 



The comparison of our results for both cases with 
the results obtained in Ref. [l3| are presented in 
Tab.lIVj 



VI. CONCLUSIONS 

We summarize our findings. Using the improved 
lattice action ([T]) and the adequate gauge fixing pro- 
cedure we completed careful study of the properties 
of the thermal color-magnetic monopoles. 

Comparing our results for thermal monopole den- 
sity and parameters of the monopole interactions with 
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b 




x' 


3.653(6) 
3.80(1) 
3.66(7) 


0.718(2) 
0.86(2) 
0.78(2) 


0.2 

0.58 
0.26 



TABLE IV: The comparison of the fit parameters for three 
cases. The Wilson case (first line) [l3|, the Symanzik case 
when loops are untouched(second line) and the Symanzik 
case when all loops eliminated(third line). 



results of Ref. [T3] we find rather small deviations, 
at the level of 10% for the density. We have found 
that this difference decreases even further when small 
dipoles are not counted. This implies universality 
for infrared thermal monopoles determined in the 
MAG. Establishing of the universality of the thermal 
monopole properties is important to prove that these 
monopoles determined after the MAG fixing are fluc- 
tuations of the gauge field relevant for infrared physics 
rather than artifacts of the gauge fixing. 

To study cutoff effects we made computations with 
two lattice spacings (using lattices with Lt = A and 
6) at temperatures T = 1.5Tc and 2Tc. We find that 
results for the thermal monopole density are not de- 
pending on the lattice spacing (see FIG. [2] and Ta- 
ble |V]) and thus they are computed in the continuum 
limit. Results for parameters of the monopole inter- 
action show slight dependence on the lattice spacing 
(see Table |lT|) and thus they are close to the contin- 
uum limit. The exception is the monopole mass which 
shows strong dependence on the lattice spacing. 

We confirmed observation made before in Ref. ^10] 
that without proper gauge fixing the systematic effects 
due to Gribov copy effects are large, e.g., up to 30% 
for the monopole density. We found that these effects 
are even more essential for the monopole trajectories 
with multiple wrappings. 

Studying the correlation functions we obtained the 
values for the Coulomb coupling constant am which 
are substantially smaller than values obtained in 
Ref. 7]. Since our results are close to the contin- 
uum limit and the same is true for results of Ref. Q 
this difference is to be explained by Gribov copy ef- 
fects. Furthermore, we found highly nonmonotonic 
behavior for both a™ and screening mass ttid near the 
transition temperature Tc. Both parameters first de- 
crease with increasing temperature and then jump up 
at the temperature 1.2Tc before monotonous depen- 



dence on T is settled. We shall admit that although 
our measurements in this range of temperature were 
done on lattices with Lt — 6, the finite lattice spacing 
effects should be checked by simulations on lattices 
with larger value of Lt to check this effect. 

Comparatively small values of am give rise to small 
values of the plasma parameter ri\i. We find that 
this parameter flattens at the value about 2 at high 
temperatures in contrast to value of 5 found in Ref. [3] . 
Still we confirm that the monopoles are in a liquid 
state at all temperatures considered. 

We have repeated the study of the monopole tra- 
jectories with multiple wrapping undertaken in Q. 
Although we obtained the values for the densities of 
such trajectories quite different from the values found 
in 0, our values for the effective chemical potential 
fj, are quite close to results of Re. 0. Moreover we 
confirm that fi goes to zero, indicating Bose-Einstein 
condensation, at the temperature very close to Tc- 

We made first studies of the effects of UV fluctua- 
tions on parameters of the monopole potential and on 
monopole thermal mass removing contributions from 
the closed loops attached to the wrapped monopole 
loops. Due to removing such contribution we were 
able to identify the currents jo of the wrapped loop 
unambiguously and thus to use all time-slices of the 
lattice in the computation of the correlators, thus de- 
creasing the statistical error substantially. 
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Appendix 

In this appendix we present a table of the densities 
for all studied temperatures. 
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